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An  exact  and  finite -element  analysis  is  used  to  predict  the  static 
bending  and  natural  vibration  response  oh  rectangular  plates  laminated 
0(i  composite  material  having  difificrent  elastic  properties  depending  upon 
whether  the  ^iber- direction  strains  are  tensile  or  compressive .  The 
analysis  is  based  on  a  theory  that  accounts  fa>r  anisotropy ,  bending - 
Stretching  coupling,  thickness  shear  deformation,  and  both  coupling 
inertia  and  rotatory  inertia.  The  finite- element  results  are  hound  to 
agree  very  closely  w.ith  the.  exact  closed- harm  solutions  in  the  case  oh  a 
cross-ply  rectangular  plate  having  simply  Supported  edges  and  subjected 
to  sinusoidally  distributed  normal- pressure  loading  or  temperature 
loading . 


Introduction 

Certain  fiber-reinforced  materials,  especially  those  with  very  soft 
matrix  material  (e.q.,  cord-rubber  composites),  exhibit,  different  elastic 
behavior  dependinq  upon  whether  the  fiber-direction  strain  is  tensile  or 
compressive  [1-31.  In  other  words,  the  tanqent.  modulus  in  tension  is 
quite  different  from  the  tanqent  modulus  in  compression.  As  a  first, 
approximation  the  uniaxial  stress-strain  behavior  of  such  materials  is 
often  represented  as  beinq  linear  with  different  slopes  (i.e.,  elastic 
moduli)  dependinq  upon  the  siqn  of  the  fiber-direction  strain.  Such  a 
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material  is  called  a  bimodulus  composite  material.  Several  macroscopic 
material  models  appropriate  for  bimodulus  fiber-reinforced  composites 
have  been  proposed  and  are  reviewed  in  [4],  It  has  been  shown  that  the 
f iber-qoverned  symmetric-compliance  model  proposed  in  [5]  aqrees  well 
with  experimental  data  for  several  materials  with  drastically  different 
elastic  properties  in  tension  and  compression. 

Prior  to  the  works  of  the  present  authors,  the  literature  available 
in  Enqlish  on  bendinq  and  vibration  analyses  of  bimodulus  plates  was 
quite  sparse  and  larqely  concerned  with  the  bending  of  bimodulus 
isotropic-material  plates.  Shapiro  [6]  considered  the  simple  problem  of 
a  circular  plate  under  a  pure  bending  moment  at  its  edge.  Kamiya  [7] 
analyzed  the  large-deflection  behavior  of  clamped  circular  plates  using  a 
finite-difference  technique.  In  [8],  Kamiya  used  the  Galerkin  method  to 
analyze  the  large-deflection  behavior  of  rectanqular  plates  under 
sinusoidally  distributed  loadinq.  The  effect  of  thickness  shear 
deformation  was  included  in  the  simple  case  of  cylindrical  bendinq  by 
Kamiya  [9].  The  only  previous  analyses  applicable  to  anisotropic 
bimodulus  material  are  the  works  of  Jones  and  Morgan  [10],  who  treated 
cylindrical  bending  of  a  thin,  cross-ply  laminate,  and  the  closed-form 
solutions  of  Kincannon  and  Bert  [11,12]  for  clamped  elliptic  plates.  All 
of  the  previous  thermoelastic  analyses  of  bimodulus-material  plates  were 
also  limited  to  isotropic  bimodulus  materials  and  midplane-symmetric 
temperature  changes  [13-17].  Apparently,  there  were  no  published  papers 
on  the  vibration  of  bimodulus  plates. 

In  the  last  two  years,  the  present  authors  have  contributed  a  series 
of  papers  to  the  literature  on  bimodulus  rectanqular  plates  [18-23]. 

These  works  are  apparently  more  general  than  any  analyses  that  have 


appeared  in  the  open  literature.  More  specifically,  these  works  differ 
from  previous  works  in  the  following  respects: 

1.  The  material  of  each  layer  is  both  elastically  and 
thermoel astical ly  orthotropic  and  bimodular. 

2.  Both  single-layer  orthotropic  and  two-layer  cross-ply  laminated 
plate  constructions  are  considered. 

3.  Thickness  shear  deformations  are  included. 

4.  Temperature  changes  through  the  thickness  as  well  as  in  the 
plane  are  considered. 

5.  Both  finite-element  and  exact  closed-form  solutions  are 

\ 

presented. 

The  objective  of  the  present  paper  is  two-fold:  to  review  the 
recent  developments  in  the  analyses  of  bimodulus-material  plates  with 
finite  dimensions  and  to  present  new  finite-element  results  for 
rectangular  plates  under  conditions  that  do  not  admit  closed-form 
solut ions . 


2.  Governing  Equations  for  Bimodulus  Plates 


In  the  following  the  equations  governing  the  shear  deformable  theory 
of  layered  composite  plates  (see  Whitney  and  Paqano  [24])  is  reviewed  and 
the  thermoelastic  constitutive  equations  of  bimodulus-material  plates  are 
presented.  In  deriving  the  qoverninq  equations  it  is  assumed  that  all  of 
the  layers  in  the  plate  remain  elastic  during  the  deformation,  the 
generalized  Hooke's  law  is  valid,  and  no  slip  occurs  between  any  two 


layers. 


Consider  a  plate  constructed  of  a  finite  number  of  uniform 


thickness,  orthotropic,  bimodulus  layers  with  arbitrary  orientations 
(i.e.,  the  material  symmetry  axes  of  each  layer  do  not  coincide,  in 
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general  with  the  plate  axes).  The  x-y  plane  lies  in  the  middle  plane,  R, 
of  the  plate  and  the  z-axis  is  normal  to  R.  The  displacement  field  in 
the  shear  deformable  theory  of  plates  is  assumed  to  be  given  by 
u(x,y,z,t)  =  uQ(x,y,t)  +  z<|»x(x,y,t) 

v(x,y,z,t)  =  vQ(x,y,t)  +  z  <py(x,y,  t )  (2.1) 

w  =  w(x,y,t) 

where  u,  v,  and  w  are  the  displacements  along  x,  y  and  z  directions, 
respectively,  u0  and  v0  are  the  in-plane  displacements  of  the  middle 
plane,  4>x  and  <|>y  are  the  bending  slopes,  and  t  is  time.  Using  the 
smal 1 -deflection  theory,  the  strain-displacement  relations  can  he 
expressed  in  the  form 


£1  =  ev  =  Ur>  v  +  Z^v  v  »  ZJ  -  ew  =  Vr,  w  +  2 ’i'w  „» 

4  a  u,a  a, a  4  y  u,y  y*y 


ec  =  y  =  u  +  v  +  z(  ii>  +  d;  ) 

8  jyw  n  \j  v  “v  v/  ^v/  v' 


'xy  o,y  o,x 


x,y  y,  x' 


Yxz  =  *x  +  w’x  ’  Yyz  =  ^y  +  w’y 
where  w,x  i  3w/3x,  etc. 

•  The  equations  of  motion  are  given  by, 
\x  +  Nfa,y  =  Puo,tt  +  R*x,tt 

N6,X  +  N2,y  =  Pv0,tt  +  R^y,tt 


e3  =  0 


(2.2) 


*1,X  +  Q2,y  =  q  +  Pw,tt 

\x  +  \y  -  Q1  =  !*x,tt  +  R  Uo,tt 

Mfc,X  +  M2,y  -  Q2  =  !*y,tt  +  R  V0,tt 

where  P,  R,  and  I  are  the  normal,  coupled-normal -rotatory,  and  rotatory 


inertia  coefficients. 


(P,R,I)  =  I 


/•  m+1 

/  (!.«.»* 


Ipw  <tz 


(?.4) 


*  z 
m  m 


p^being  the  material  density  of  the  m-th  layer,  q  is  the  transversely 


distributed  force,  and  N. ,  Q-,  and  M-  are  the  respective  inplane  and 

transverse  stress  and  moment  resultants  defined  by 

h/2  h/2 

(N-.M  )  =J  (l,z)o.  dz  (i  =  1,2,6)  ,  (QpQ2)  =J  (o5,o4)dz, 

(2.5) 

where  h  is  the  plate  (laminate)  thickness,  and  the  so-called  contracted 
subscript  notation  is  employed  to  denote  the  stress  components. 

Assuminq  that  the  only  plane  of  symmetry  existing  is  in  the  plane  of 
the  plate,  the  thermoelastic  constitutive  relations  for  each  layer  (t) 
are  taken  to  be  monoclinic  and  bimodular  as  follows: 


'°nia  Qi2k£  0 

j  "I  -  <W 

°2 

= 

Ql2|a  Q22k£  ° 

1  -  a2k/ 

°b 

-  0  0  Qfefckt- 

(  Eb 

{:;!■[ 


c<*5ia 
C4bia  C5Sk* 


]|:;1 


(2.7) 


Here  T  is  the  temperature  change  measured  from  the  strain-free 
temperature,  the  C's  are  Cauchy  elastic  shear  stiffnesses,  the  Q's  are 
plane-stress  reduced  stiffnesses,  and  the  a's  are  thermal-expansion 
coefficients.  The  third  subscript  in  Q..^  and  C-^  (and  second  in 
ct -ke)  refers  to  the  bimodular  characteristics:  k  =  1  denotes  properties 
associated  with  fiber-direction  tension,  k  =  2  denotes  fiber-direction 
compression.  Also,  subscript  l  refers  to  the  individual  layer  number, 
i.e.,  i  -  1  and  2  for  a  two- layer  laminate. 

Substituting  (2.2),  (2.6)  and  (2.7)  into  (2.5),  we  obtain  the 


constitutive  equations  for  an  arbitrarily  laminated  plate. 
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[A] 

JB1 

[bit 

[0] 

M- 

({«]) 

(§0  •  [is  *s]  {y  • 


P  ■  ISI 


where 


o 

e,  =  u 
l  o,x 

<  -  =  ili 

1  *x,x 


°  =  v 

2  o,y 


*9  = 


y,.y 


cl  =  uo,y  +  vo,x 
k  x,y  rv,x 


(2.8) 

(2.9) 

(2.10) 


and  A..,  B--,  D-.,  S,,  are  the  respective  inplane,  bendinq-inplane 

^  J  '  J  J  ^  J 

coupling,  bending  or  twistinq,  and  thickness-shear  stiffnesses  defined  as 
f o 1 1 ows : 

,h/2 

^Aij’Bij’Dij^  =/_h/,2^1,z’zi>^ijkJt  dz  =  1>2’6) 


(2.11) 


h/2 


Sij  ‘  K< 


iKj/_h//?  ^ijk*  dz 


The  stress  and  moment  resultants,  nT  and  mT,  due  to  thermal  loadinq,  are 
defined  by 

(nI.mI)  =  zflm+]  £  Qjjkt  «jlu(VzTi)dz  ,  (ij  =  1,2,6)  (2.12) 

where  the  temperature  change  from  a  reference  state  is  assumed  to  be 
linear  through  the  thickness,  consistent  with  the  plate  theory: 

T(x,t,z)  =  TQ(x,y)  +  zTj(x,y).  (2.13) 

This  completes  the  description  of  the  field  equations.  Closed-form 
solutions  to  these  equations  are  possible  in  the  particular  case  of  a 
cross-ply  plate  whose  stiffness  coefficients  (A.^,  B.j  and  0.^)  with 
subscripts  16,  26,  and  45  are  zero,  and  whose  edges  are  hinged  flexural ly 
but  free  to  move  in  a  direction  normal  to  each  edqe.  The  temperature 
distribution  and  the  normal -pressure  loadings  must  be  sinusoidal  in  order 
to  obtain  the  closed-form  solution  (see  [18-23]). 


previous  section  follows  along  the  same  lines  as  for  ordinary  (i.e.,  not 

bimodulus)  plates  (see  Reddy  [25]).  The  main  difference  lies  in  the 

computation  of  the  stiffness  coefficients,  which  depend  on  the  stresses 

in  the  plate.  We  summarize  the  results  in  the  following. 

Consider  a  f inite-element  representation  of  the  midplane  of  the 

plate  by  quadr i lateral  isoparametric  elements.  Over  a  typical  element, 

Re,  the  displacements  (u  ,  v  ,w,i|i  ,  ip  )  are  approximated  by  expressions 

oo  x  y 


uo  =  f,  Vi  •  vo  =  ?  Vi’  etc- 


(3.1) 


where  <t>.(x,y)  (i  =  l,2,...,n)  are  interpolation  functions  associated  with 
the  linear  (n  =  4)  or  quadratic  (n  =  8  or  9)  rectangular  element. 
Substituting  (3.1)  into  the  weak  form  associated  with  (2.3),  (2.8),  and 


(2.9),  we  obtain 

([«]  -  o»2[M]){a}  =  {F) 


(3.2) 


The  elements  of  the  stiffness  matrix  [K]  and  mass  matrix  [M]  are  given  by 


■  AuG!j  *  \6Gij 

K‘j  ■  BnGu  ♦  B«Glj 
^  ^ 
KtJ  '  B«SU  *  B22G)j 


K|j2  *  a,2=1]  *  • 

I'D  *  8.^  *  • 

Ki.i  *  *  8 12®j?  • 

"D  ■  ^  • 

‘8  •  ^  ■ 


KU  ■  Dll«?j  *  °66GU  *  SS5Gij  •  I'D  *  0 , 2®  i  j  *  »t*G)i  • 

I'D  *  *  S,,^  ,  K{j  *  ■  o  . 


(3.3) 
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ij 

ij 
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R 

Sn  _ 
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ij 

(3.4) 

(3.5) 

(3.6) 


f,  =  n;  +  n'  ,  f,  =  »;  +  n,  ,  f,  =  q  , 

fu  =  mI  +  m!  ,  f,  =  mI  +  .  (3.7) 

4  l,x  fe,y  5  6,x  2,y  ' 

In  the  present  study  the  4-node  and  8-node  rectangular  isoparametr ic 
elements  are  employed. 


4.  Numerical  Results  and  Discussion 

Numerical  results  for  deflections,  neutral -surface  locations,  and 
vibration  frequencies  are  presented  for  plates  under  various  edqe 
conditions  and  loadinqs.  The  discussion  of  the  results  is  divided  into 
three  parts  according  to  bending  due  to  normal  pressure,  bendinq  due  to 
thermal  loading,  and  free  vibration.  In  all  three  parts,  numerical 
results  are  presented  first  for  plate  problems  which  admit  closed-form 
solutions.  All  of  the  finite-element  results  presented  herein  were 
obtained  on  an  IBM  370/158  computer  using  the  double-precision 
arithmetic. 


Table  1 .  Elastic  properties  for  two  ti record-rubber,  unidirect iona  1_, 

bimodulus  composite  materials 


Aram id  1 

Rubber  (AR) 

Polyester- 

-Rubber  (PR) 

/ 

Property  and  Units 

Tension 

(k=l) 

Compression 

(k=2) 

Tension 

(k=l) 

Compression 

(k=2) 

Lonqitudinal  Young's 

modulus,  GPa  (E1}) 
Transverse  Younq's  modulus. 

3.58 

0.0120 

0.617 

0.0369 

GPa  (E  2) 

Major  Poisson's  ratio. 

0.00909 

0.0120 

0.00800 

0.0106 

dimensionless3  (vi;) 
Longitudinal-transverse 

0.416 

0.205 

0.475 

0. 185 

shear  modulus,  GPab  (G12) 
Transverse-thickness  shear 

0.00370 

0.00370 

0.00267 

0.00267 

modulus,  GPa  (G2j) 

0.00290 

0.00499 

0.00233 

0.00475 

aThe  minor  Poisson's  ratio  is  assumed  to  be  given  by  the  reciprocal 
relation. 

bThe  lonqitudinal-thickness  shear  modulus  is  assumed  to  be  equal  to  the 
longitudinal-transverse  shear  modulus. 


Static  Bending  Due  to  Normal  Pressure 

To  show  the  accuracy  of  the  finite-element  solutions,  a  comparison 
is  made  of  the  present  results  with  the  closed-form  results  obtained  in 
[21].  Table  2  contains  the  neutral-surface  locations,  and  the 
nondiinensional ized  transverse  deflection  of  single-layer  and  two-layer 
(0°/90'’),  simply-supported  (SS),  rectangular  plates  ( arainid-rubber )  under 
sinusoidally  (SSL)  distributed  transverse  load,  q  =  qQsin( nx/a)sin( ny/b) . 
As  can  be  seen  from  Table  2,  the  agreement  between  the  closed-form 
solution  (CFS)  and  finite-element  solution  (FES)  is  extremely  good.  It 
should  be  noted  that  if  the  material  is  treated  as  an  ordinary  material, 
with  material  properties  taken  as  the  averaqe  of  compressive  and  tensile 
properties,  the  deflections  will  be  under-estimated  siqnif icantly. 


Table  2.  Neutral-surface  locations  and  dimensionless  center  deflections  for  simply-supported 
thick  rectangular  plates  of  aramid-rubber  (b/h  =  10,  shear  correction  factor  =  5/6) 
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The  aramid-rubber  plates,  in  both  the  single-  and  cross-ply  cases, 
have  noticeably  larger  values  of  Zx  than  the  polyester-rubber  ones,  due 
to  the  more  pronounced  bimodulus  effect  in  the  fiber-direction  Young's 
modulus  of  the  aramid-rubber.  There  are  only  very  slight  differences  in 
Zx  and  deflection  in  going  from  a  single-layer  plate  to  a  cross-ply 
laminate.  This  is  in  considerable  contrast  to  ordinary  materials. 

Having  validated  the  finite-element  model,  the  model  is  then  used 
for  other  combinations  of  loading  and  boundary  conditions  that  are  not 
amenable  to  closed-form  solutions.  Figure  1  shows  plots  of  the 
nondimensional ized  deflection  versus  the  side-to-thickness  ratio  for  two- 
layer  ( 0*790° )  simply-supported  and  clamped  (CC)  square  plates  under 
various  loadings.  The  effect  of  the  thickness  shear  on  the  deflections 
is  apparent  from  the  graphs. 

Static  Bending  Due  to  Thermal  Loading 

Next,  numerical  results  are  presented  for  static  bending  of  plates 
under  thermal  gradients  only  (i.e.,  Tq  =  0).  The  same  orthotropic 
bimodulus  materials  as  listed  in  Table  1  are  used.  Since  no  measured 
values  of  thermal-expansion  coefficients  of  these  materials  are 
available,  the  following  ratios  of  the  thermal -expansion  coefficients  are 
used: 

a^/otj  =  0.5  ;  <*2/012  =1.0  ;  aVaz  =  0.1.  (4.3) 

To  show  the  accuracy  of  the  finite-element  solutions,  a  comparison 
is  made  of  the  present  results  with  the  closed-form  results  obtained  in 
r?2].  Table  3  contains  the  neutral-surface  locations,  and  the 
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nondimensional ized  transverse  deflection  of  sinqle-1 ayer  and  two-layer 
(0°/90°),  simply-supported  (SS),  rectangular  plates  (polyester-rubber) 
under  sinusoidally  (SST)  distributed  temperature,  T  =  T  xa  sin(nx/a) 
sin(wy/b).  As  can  be  seen  from  Table  3,  the  agreement  between  the 
closed-form  solution  (CFS)  and  finite-element  solution  (FES)  is  extremely 
good  (see  also  Fiq.  2) . 

Table  3.  Neutral-surface  locations  and  dimensionless  center  deflections 
for  simply-supported  thick  (b/h  =  10,  Kj  =  K2  =  5/6) 
rectangular  plates  of  polyester -rubber  subjected  to  sinusoidal 
temperature  gradient,  T  =  zT  l  sin(  iTx/a)sin(  ny/b) . 


a/b 

single- 

layer 

two- layer  + 

Zx  (" 

zx/h) 

w 

f+ 

Zx  (  = 

zx/h) 

w 

CFS 

FES 

CFS 

FES 

CFS 

FES 

CFS 

FES 

0.50 

0.1031 

0.1030 

0.482 

0.482 

0.2599 

0.2541 

0.928 

0.894 

0.75 

0.1184 

0.1183 

1.157 

0.158 

0.2554 

0.2436 

0.199 

0.187 

1.00 

0.1308 

0.1308 

2.160 

2.161 

0.2398 

0.2294 

0.309 

0.296 

1.25 

0.1360 

0.1360 

3.410 

3.409 

0.2119 

0.2035 

0.382 

0.382 

1.50 

0.1332 

0.1331 

4.737 

4.736 

0.1734 

0.1679 

0.443 

0.436 

1.75 

0.1234 

0.1233 

5.975 

5.970 

0.1284 

0.1252 

0.472 

0.466 

2.0 

0.1078 

0.1078 

7.024 

7.017 

0.0815 

0.0804 

0.489 

0.483 

+Neutral -surf ace  location  Zy  is  not  listed  here  (see  [22]) 

++w  =  (wh/a*  Tj  b2)  10 

Next,  finite-element  results  are  presented  for  other  combinations  of 
loading  and  boundary  conditions  that  are  not  amenable  to  closed-form 
solutions.  Figure  3  shows  plots  of  the  nondimensional ized  deflection 
versus  the  side-to-thickness  ratio  for  two-layer  (0'/90°)  and  one-layer 
simply-supported  (SS)  square  plates  (polyester-rubber)  under  uniform 
temperature  gradient  through  the  thickness.  The  effect  of  the  thickness 
shear  on  the  deflections  is  apparent  from  the  graphs. 


14 


Fig. 3  Nondimens i onal i zed  center  deflection  vs.  side-to-thickness 
ratio  for  two-layer  simply  supported  (SS)  square  plates 
under  uniformly  distributed  temperature  gradient  through 
the  thickness. 


Natural  Vibration  Results 

Before  we  present  the  numerical  results  of  free  vibration,  we 
discuss  the  effect  of  bimodulus  action  on  plate  stiffness  in  different 
portions  of  a  vibration  cycle.  First  consider  a  single-layer,  bimodulus 
material  plate  at  the  two  extremes  of  its  deflection  (Fiq.  4).  During 
the  first  half  cycle  (see  Fig.  4a)  the  top  surface  is  in  compression  and 
the  bottom  in  tension,  thus  causinq  the  neutral  surface  for  ex  to  be 
positive  (znx  >  0),  i.e.,  a  certain  distance  below  the  plate  midplane. 
During  the  second  half  cycle  (Fiq.  4b)  the  top  surface  is  in  tension  and 
the  bottom  is  in  compression,  just  opposite  to  that  encountered  in  the 
first  half  cycle.  However,  the  magnitude  of  znx  is  the  same  as  in  the 
first  half  cycle.  Thus,  it  can  be  concluded  that  the  effective  stiffness 
(and  thus  the  frequency)  associated  with  the  second  half  cycle  is 
identical  to  that  of  the  first  half  cycle  and  either  modal  shape,  Fiq.  4a 
or  4b,  will  give  the  same  computational  result  for  the  natural 
frequencies . 

Now  consider  a  two-ply  laminate  with  the  bottom  layer  (layer  l  =  1) 
oriented  at  0°  and  the  top  layer  (l  =  2)  at  90’  (see  Fiq.  4).  Initially, 
as  shown  in  Fig.  4c,  the  neutral  surface  for  ex  falls  below  the 
interface,  within  the  O'  layer,  while  the  neutral  surface  for  falls 
above  the  interface,  completely  within  the  90’  layer.  In  the  latter 
portion  of  the  cycle.  Fig.  4d,  the  ex  neutral  surface  falls  outside  of 
the  layer,  and  the  neutral  surface  falls  outside  of  the  90’  layer. 

Therefore,  compressive  properties  are  used  for  the  entire  0°  layer,  and 
tensile  ones  for  the  90’  layer.  Thus,  for  a  two-layer  cross-ply  laminate 
the  plate  stiffnesses  actinq  in  the  two  portions  of  a  cycle  are  different 
and  the  associated  frequencies  are  also  different  (except  for  a  square 
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(A)  FIRST  HALF  CYCLE 


(B)  SECOND  HALF  CYCLE 


SINGLE-LAYER,  0°  BIMODULUS  PLATE 


(C)  FIRST  PORTION  (D)  SECOND  PORTION 

OF  CYCLE  OF  CYCLE 

TWO-LAYER,  CROSS-PLY  PLATE 
(BOTTOM  LAYER  AT  0°,  TOP  LAYER  AT  90°) 


Fig.  4  Bimodulus  action  in  single-layer  and  cross-ply  plates  in 
the  fundamental  mode  of  vibration.  Cross-hatched  regions 
denote  tension  in  the  respective  fiber  directions. 
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plate).  Denoting  the  frequencies  associated  with  the  two  portions  of  a 
vibration  cycle  by  Wj  and  ui2,  we  can  easily  show  that  the  averaqe 
frequency  (u>)  over  the  entire  cycle  is 

of1  =  ( 1  /  2 )  ( w  j 1  +  o'1)  (4.4) 

Thus,  the  computational  procedure  required  for  a  cross-ply  plate  is  to 
calculate  Wj  and  u>2  associated  with  modal  shapes  shown  in  Figs.  4c  and 
4d,  respectively,  and  then  to  apply  equation  (4.4). 

The  question  arises  as  to  the  possibility  of  a  discontinuity  in 
energy  at  the  instant  of  time  between  the  two  portions  of  a  cycle. 
However,  it  was  shown  in  [23]  that  there  is  no  such  discontinuity. 

Numerical  results  are  presented  for  relatively  thick  plates  (b/h  = 
10)  of  aramid-rubber.  Plate  aspect  ratio  affects  the  fundamental 
frequency  much  more  than  it  does  the  neutral-surface  portion  (see  Tables 
4-5).  For  the  cross-ply  case,  only  and  ^  is  close  to  Z^.  In 

XX  y 

Fig.  5  the  good  agreement  between  the  CFS  and  FES  results  is  clearly 
shown. 

Since  there  are  very  drastic  changes  in  neutral-surface  locations 
(for  example,  znx  for  aramid-rubber  goes  from  approximately  0.4h  to  - 
0.03h)  from  one  cycle  portion  to  the  other,  there  is  a  question  that  may 
arise  regarding  a  possible  transient  action.  However,  the  neutral 
surfaces  are  just  boundaries  between  the  tensile  and  compressive  reqions 
(analogous  to  the  elastic-plastic  boundary  in  elastoplastic  problems). 
Thus,  they  have  no  mass  and  there  is  no  transient  action. 


Table  4  Dimensionless  fiber-direction  neutral-surface  locations  for  single-layer  (0°)  and  two 
layer  (0*/90*)  plates  of  polyester-rubber  material  (b/h  =  10). 
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Table  5  Dimensionless  fundamental  frequencies  in  the  first  partial 

cycle,  second  partial  cycle,  and  complete  cycle  of  motion  for 
two-layer,  cross-ply  plates  having  b/h  =  10  by  closed  form  and 
finite  element  methods. 


a/b 

u>1b2(P/E^2h3)1/2 

u>^b2(P/E^2h3) 1/2 

ii)2(P/E^2t 

,3)172 

CFS 

FES 

CFS 

FES 

CFS 

FES 

Aramid-Rubber: 

0.5 

19.38 

20.23 

13.88 

14.55 

16.18 

16.93 

0.7 

11.60 

12.17 

9.353 

9.807 

10.35 

10.86 

1.0 

7.038 

7.386 

7.038 

7.364 

7.038 

7.375 

1.4 

4.838 

5.045 

6.037 

6.356 

5.371 

5.625 

2.0 

3.712 

3.909 

5.551 

5.821 

4.449 

4.677 

Polyester-Rubber : 

0.5 

19.12 

19.81 

15.95 

16.61 

17.39 

18.07 

0.7 

11.43 

11.92 

10.04 

10.45 

10.69 

11.14 

1.0 

7.084 

7.406 

7.085 

7.394 

7.085 

7.400 

1.4 

5.164 

5.407 

5.928 

6.193 

5.520 

5.773 

2.0 

4.310 

4.518 

5.435 

5.688 

4.807 

5.036 

5.  Summary  and  Conclusions 

Both  closed-form  and  finite-element  solutions  are  presented  for 
static  bending  and  free  vibrations  of  single-layer  and  two-layer  cross- 
ply  plates  of  bimodulus  composite  materials.  In  the  case  of  simply 
supported  plates  under  sinusoidal  temperature  distribution,  normal- 
pressure  loading,  and  free  vibration,  the  closed-form  and  the  finite- 
element  solutions  for  neutral-surface  locations,  center  deflections,  and 
fundamental  frequencies  are  in  excellent  agreement.  The  finite-element 
method  is  then  used  to  obtain  solutions  for  simply-supported  plates  under 
uniform  loading,  for  which  no  closed-form  solution  can  be  obtained  for 
the  bimodulus  case. 

The  research  reported  here  has  recently  been  extended  to 
mechanically  loaded  thin  cylindrical  shells  [26],  free  vibration  of  thin 
and  thick  cylindrical  shells  [27],  and  thermally  and  mechanically  loaded 
thick  cylindrical  shells  [28]. 
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theory  that  accounts  for  anisotropy,  bending-stretching  coupling,  thickness 
shear  deformation,  and  both  coupling  inertia  and  rotatory  inertia.  The  finite 
element  results  are  found  to  agree  very  closely  with  the  exact  closed-form 
solutions  in  the  case  of  a  cross-ply  rectangular  plate  having  simply  (over^ 
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IB.  Supplementary  Notes  -  Cont'd 

of  Layered  Composite  Plates  and  Shells  of  Bimodulus  Materials”  to  be 
presented  at  the  Second  International  Conference  on  Numerical  Methods  in 
Thermal  Problems,  Venice,  Italy,  July  7-10,  1981,  and  will  appear  in  its 
Proceedings. 
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supported  edges  and  subjected  to  sinusoidally  distributed  normal -pressure 
loading  or  temperature  loading. 
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